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ABSTRACT 



Using galaxy clusters from the ESO Distant Cluster Survey, we study how the dis- 
tribution of galaxies along the colour-magnitude relation has evolved since z ~ 0.8. 
While red-sequence galaxies in all these clusters are well described by an old, pas- 
sively evolving population, we confirm our previous finding of a significant evolution 
in their luminosity distribution as a function of redshift. When compared to galaxy 
clusters in the local Universe, the high redshift EDisCS clusters exhibit a significant 
deficit of faint red galaxies. Combining clusters in three different redshift bins, and 
defining as 'faint' all galaxies in the range 0.4 >L/L* > 0.1, we find a clear decrease 
in the luminous-to-faint ratio of red galaxies from z ~ 0.8 to z ~ 0.4. The amount 
of such a decrease appears to be in qualitative agreement with predictions of a model 
where the blue bright galaxies that populate the colour-magnitude diagram of high 
redshift clusters, have their star formation suppressed by the hostile cluster environ- 
ment. Although model results need to be interpreted with caution, our findings clearly 
indicate that the red-sequence population of high redshift clusters docs not contain 
all progenitors of nearby red-sequence cluster galaxies. A significant fraction of these 
must have moved onto the red-sequence below z ~ 0.8. 



Key words: galaxies: clusters: general - galaxies: evolution - galaxies: fundamental 
parameters - galaxies: luminosity function, mass function. 



1 INTRODUCTION 



similar environments over a relatively long time baseline. 
In addition, they offer the practical advantage of provid- 
ing many galaxies in a relatively small region of the sky 
and all approximately at the same redshift. This allows effi- 
cient observation even with modest fields of view and mod- 
est amounts of telescope time. It should be noted, however, 
that in order to establish that physical processes related to 
the cluster environment are indeed playing a role, it is nec- 
essary to compare the evolution of similar galaxies in dif- 
ferent environments (i.e. in the clusters and in the 'field'). 



Galaxy clusters may be considered as laboratories for study- 
ing the physical processes that drive galaxy evolution. They 
offer the possibility to trace the properties of galaxies in 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, as part of large programme 166.A-0162 (the ESO 
Distant Cluster Survey) 
f Email: gdelucia@mpa-garching.mpg.de 
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In addition, galaxy clusters represent a biased environment 
for evolutionary studies. In the current standard cosmogony, 
clusters originate from the gravitational collapse of the high- 
est (and rarest) peaks of primordial density perturbations, 
and evolutionary processes in these regions occur at an ac- 
celerated pace with respect to regions of the Universe with 
average density. 

The technical capabilities achieved in recent years 
have provided a rapidly growing dat a base on high- 



200 ll IValtchanov et alJ 12004 


iGladders & Yed 2005: 


Kodama et alJ 


2005; 


IWhite et alJ 


2005; 


IWilson et al.ll2006). 



These, interpreted using the latest theoretical techniques 
JCole et al.ll2000l: lHatton et al-lliool ; ISpringel et all 120051: 

iDe Lucia et al.ll2006l) . should provide important constraints 
on the physical mechanisms driving the formation and the 
evolution of cluster galaxies. 

Early studies of the galaxy population in distant 
clusters pointed o ut significant differences with respect 
to nearby systems feutcher fc Oemlerlll98 4) . More recent 
work has provided us with a much more detailed pic- 
ture of these di fferences (a very incomplete list of recent 
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2005; 


IWhite et alJ I2005 1 


Poggianti et al. 


20061: IStrazzuilo et alJ 


2006). One interesting outcome of 



these studies has been the discovery that a tight relation 
between the colours and the magnitudes of bright ellipti- 
cal galaxi es holds up to the highest redshifts probed so 
far iBlakeslee et alJ20ollDe Lucia et all2004l ; lrlolden et alJ 
l2004l;lMei et al.ll2006ft . 

The existence of a colour-magnitude relation 
(hereinafter CMR) has been known for a long time 
ide Vaucouleursl [l96ll IVisvanathan fc Sandagel Il977l) . At 
least in nearby clusters dde Propris et al.lll998l) . it appears 
to extend 5 — 6 mag faint-ward of the Brightest Cluster 
Galaxy (BCG). At the present epoch, the CMR can 
be interpreted either as a result of differing ages (bluer 
galaxies being younger), or of differing metallicities (bluer 
galaxies being more metal poor), or a s a combination of the 
two (Ferreras, Chario t & Silk 1999; iTerlevich et~al"lll999l 
IPoggianti et al.l|2001al . and references therein). 

The mere existence of a tight relation at high red- 
shift favours the metallicity interpretation and naively ap- 
pears to make the age explanation untenable. The reason 
for this is that if one assumes that all present-day red- 
sequence galaxies are still identified as red-sequence mem- 
bers in high-redshift clusters, then if the CMR were primar- 
ily age driven it would change dramatically with increasing 
redshift as small galaxies approach their fo rmation epoch 
and p rogressively become brighter and bluer ( Kodama et al. 
1998). This expectation is in contrast with observational re- 
sults which show that the slope of the CMR does not c hange 
appreciably over the redsh ift interval — 1 ( Gladdcrs et alJ 
ll998l:IStanford et al.lll99St and more recent work mentioned 
above) . 

One possible simple interpretation is then that cluster 
elliptical galaxies represent a passively evolving population 
formed at high-redshift (z > 2 — 3) in a short duration event 
(but see the discussion by iBower. Kodama fc Terlevichl 
1998). In this scenario - often referred to as monolithic - 
the CMR arises through the effects of supernovae winds: su- 
pernovae explosions heat the interstellar medium triggering 



galactic winds whenever the thermal energy of the gas ex- 
ceeds its gravitational binding energy. Since smaller galaxies 
have shallower potential wells, this results in greater mass 
loss by smaller systems, naturally establishing the observed 
CMR. A difficulty may be that the observed CMR shows 
no sign of a t urnover at h igh mass of the kind predicted by 
such models IILarsorJll974h . 

This simple model may be too naive for explaining the 
origin and the evolution of the observed CMR. In the mono- 
lithic scenario a galaxy has a single well-defined progenitor 
at each redshift and its evolution is described by simple 
smooth functions of time. This is not true in the current 
standard cosmological paradigm, where a single galaxy to- 
day corresponds to the ensemble of a ll its progenitors at 
any p revious redshift (see discussion in IDe Lucia fc Blaizod 
2006). It is not obvious then that high-redshift red-sequence 
galaxies contain all or even most of the progenitors of nearby 
red-sequence cluster galaxies. Indeed the results we present 
below give direct evidence that this is not the case. In addi- 
tion, the simple model described above, clearly neglects the 
infall of "new" galaxies during cosmological growth of the 
cluster. If the intra-cluster environment is associated with 
suppression of star formation, then these galaxies would 
become redder and fainter and might also join the red- 
sequence galaxy population at lower redshifts. 

An alternative scenario has been pro- 
posed bv iKa uffmann fc Chari ot] ( Tl 998) - see also 
IDe Lucia. Kauffrmtmi^^WIiitel i2004T) - in the frame- 
work of hierarchical models of galaxy formation. In these 
models, the CMR arises as a result of the fact that more 
massive ellipticals originate from the mergers of more 
massive - and more metal-rich - disk systems. The models 
show a well defined red-sequence, still mainly driven by 
metallicity differences, that is in place up to redshift ~ 2, 
although with a scatter that is larger than that observed. 

Recent work on the observed colour-magnitude rela- 
tion of high-redshift clusters has pointed out a new and 
still controversial result concerning an apparent 'truncation ' 
of the CMR at red shift about 0.8 JPe Lucia et alJ 120041 : 
iKodama et~ai]|200^ . 

In iPe Lucia et~aT1 (12004) . we analysed the colour- 
magnitude relation of 4 clusters in the redshift interval 
0.7-0.8 from the ESO Pistant Cluster Survey (hereinafter 
EPisCS) and found a deficiency of low luminosity passive 
red galaxies with respect to the nearby Coma cluster. A de- 
crease in the number of faint red galaxies was detected in 
all clusters under investigation but one (with low number 
of cluster members) , although the significance of the deficit 
was only at about the 3<r level. In this paper we extend 
our analysis to the full EPisCS sample covering the redshift 
range 0.4 — 0.8 and a wide range of structural properties. 
The plan of the paper is as follows. The observational data 
used for our study are briefly described in Sec. |2] In Sec. |H] 
we present the criteria used to define cluster membership 
and in Sec. [3] we present the colour-magnitude relation for 
all the clusters in the EPisCS sample. In Sec.|S]we study the 
distribution of galaxies along the red-sequence, and discuss 
its dependence on redshift and on cluster velocity dispersion. 
The red-sequence galaxy distribution in nearby clusters is 
studied in Sec.|S| In Sec.|7| we interpret the evolution mea- 
sured as a function of redshift in terms of simple population 
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synthesis models. Finally, in Sec. |HJ we discuss our results 
and give our conclusions. 

Throughout this paper we will assume a ACDM cosmol- 
ogy: H = 70kms" 1 Mpc -1 , fi m = 0.3 and fi A = 0.7. With 
this cosmology, z ~ 0.8 - the highest redshift probed by our 
cluster sample - corresponds to more that 50 per cent of the 
look-back time to the Big Bang. Throughout this paper we 
use Vega magnitudes, unless otherwise stated. 



2 THE DATA 

EDisCS is an ESO Large Programme aimed at the study of 
cluster structure and cluster galaxy evolution over a signifi- 
cant fraction of cosmic time. The complete EDisCS dataset 
provides homogeneous photometry and spectroscopy for 20 
fields containing galaxy clusters at z = 0.4-1. Clusters can- 
didates were selected f rom the Las Campana s Distant Clus- 
ter Survey (LCDCS) of lGonzalez et al.l fcOOll) by identifying 
surface brightness excesses using a very wide filter (4500- 
7500 A) in order to maximise the signal-to-noise of distant 
clusters against the sky. The EDisCS sample of 20 clusters 
was constructed selecting 30 from the highest surface bright- 
ness candidates in the LCDCS, and confirming the presence 
of an apparent cluster and of a possible red sequence with 
VLT 20min exposures in two filters. From these 30 candi- 
dates, we then followed up 10 among the highest surface 
brightness clusters in the LCDCS in each of the ranges of 
estimated redshift 0.45 < z cs t < 0.55 and 0.75 < z cst < 0.85. 
In the following, we will often refer to these as the interme- 
diate and high redshift samples respectively. 

As a consequence of the scatter of the estimated red- 
shifts about the true value, we have ended up with a set of 
clusters distributed relatively smoothly between z = 0.42 
and z — 0.96, rather than two samples concentrated at 
0.5 and 0.8, as originally planned. De tails on the s e lectio n 
of cluster candidates can be found in IWhite et ail £2005). 
Our follow-up progra mme obtained deep optical photome- 
try with FORS2/VLT dWhite et alJl2005l) . near-IR photom- 
etry with SOFI/NTT (Aragon-Salamanca et al. in prepara- 
tion), and multislit spectroscopy with FORS2/VLT for the 
20 fields (Halliday et al. 2004; Milvang- Jensen et al. in prepa- 
ration). ACS/HST mosaic imaging of 10 of the highest red- 
shift clusters has also been acquired (Desai et al. in prepa- 
ration). For three E DisCS clusters, n arrowband H a imaging 
has been acquired llFinn^^^J|200 5 ) and for three clusters 
we have XMM data \ Johnson et al]|200rJ) . 

The optical ground-based photometry and a first basic 
character isation of our samp le of clusters as a whole, is pre- 
sented in lWhite et alJ (l2005h . In brief, our optical photome- 
try consists of V, R, and I imaging for the 10 highest redshift 
cluster candidates and B, V, and I imaging for the remain- 
ing 10 intermediate redshift cluster candidates. Total inte- 
gration times were typically 45 minutes at the lower redshift 
and 2 hours at the higher redshift. Object catalogues have 
been created using the SExtractor software version 2.2.2 
jBertin fc Arnoutsl llOQe) in 'two-image' mode using the I- 
band images as detection reference images. Magnitudes and 
colours have been measured on the seeing-matched images 
(to 0'.'8 - the typical seeing in our IR images) using fixed 
circular apertures. Throughout this paper, we correct mag- 
nitudes and colours for Galactic extinction according to the 



maps of lSchlegel. Finkbeiner fc Davisl lll998l) and a standard 
Milky Way reddening curve. We refer to lWhite et al] (l2005h 
for details about our photometry. In the following, we will 
use magnitudes and colours measured using a fixed circular 
aperture with l'.'O radius. This choice has been adopted to 
simplify the comparison with the Coma cluster, as we will 
explain later in the paper. The cluster velocity dispersions 
and R200 we use in t he following are the same as used in 
iPoggiaritT et al.1 (120061) and are listed in Table 1 of that pa- 
per. 

As a part of our programme, we also obtained spectra 
for > 100 galaxies per cluster field (typical exposure times 
were 4 hours for the high redshift candidates and 2 hours for 
the intermediate redshift clusters). The spectroscopic selec- 
tion, observations, dat a reduction, and spec troscopic cata- 
logues are presented in lHallidav et all <l2004li and Milvang- 
Jensen et al. (in preparation). As explained in lWhite et al] 
(2005), deep spectroscopy was not obtained for two of the 
EDisCS fields (C11122.9-1136 1 and cll238.5-1144 2 ), which 
are not included in the present study. 



3 CLUSTER MEMBERSHIP 

Although complications arise from the existence of redshift 
space distortions, spectroscopic redshifts provide the opti- 
mal technique to determine cluster membership. However, 
obtaining spectroscopic redshifts for large numbers of faint 
objects is not feasible within the available time with cur- 
rent instrumentation, even for samples just beyond z = 0.5. 
A standard method to correct for field contamination, in 
absen ce of spectroscopy, is to use st atistical field subtrac- 
tion dAr agon-Sa lamanca et al] Il993l : IStanford et al.l Il998l : 
iKodama fc BoweJl200lR 7 a 'cluster-free' field is used to de- 
termine the number of contaminating galaxies as a function 
of magnitude and/or colour. Th is method becomes increas- 
ingly uncertain at high redshift: iDriver et al] dl998l) . for ex- 
ample, used simulations to show that it is already unreliable 
at z > 0.3. In addition, this approach does not provide the 
likelihood of being a cluster member on a galaxy-by-galaxy 
basis. 

In the last decade, the techniques used to deter- 
mine photometric redshifts have become much more pre- 
cise, suggesting that they can be used to address specific 
scientific questions ]Bemtez| |2000j: [Bolzonella et al.l l200d : 
iRudnick et aHl200ll: lFirth et al]l2003l) . An important by- 
product is an estimate of the spectral type for each observed 
galaxy. Errors in estimated photometric redshifts are much 
larger than typical errors in spectroscopic redshifts. In addi- 
tion, systematics or degeneracies are often present because 
of uncertainties in the redshift evolution of spectral energy 
distributions and/or insufficient calibrating spectroscopy for 
the magnitude range sampled by the photometric data. 

Given the difficulties mentioned above, we decided in 
the present study to use both a 'classical' statistical field 



1 Only one short exposure mask was obtained for this field, show- 
ing no evidence of a concentration of galaxies at any specific red- 
shift. 

2 Only two short exposure masks were obtained for this field for 
which we also do not have NIR data. 
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subtraction and a membership criterion based on photomet- 
ric redshift information. We give a brief description of both 
methods in the following sections. 

3.1 Photometric redshifts 

Photo metric redshifts were computed using two different 
codes jBolzonella et alJl200fi iRudnick et alJl200lt) in order 
to provide better control of the systematics in the identifica- 
tion of likely non-members. The two codes employed in this 
study are based on the use of similar SED fitting procedures, 
but different template spectra, different minimisation algo- 
rithms, and a number of other different details (see the orig- 
inal papers) . The performance of these codes on the EDisCS 
dataset will be examined in Pello et al. (in preparation). 

For the purposes of this analysis, the codes were run al- 
lowing a maximum photometric redshift of 2 and assuming 
a 5 per cent minimum flux error for the photometry. Where 
they can be checked, the photometric redshifts of the galax- 
ies in our sample are quite accurate with < |z spC c — z P hot >= 
0.06-0.08. There is no systematic trend between 2 p hot — z Bpcc 
and z S pec and the percentage of catastrophic failures, i.e. the 
fraction of objects with \z spcc — z p hot > 0.3, is of the order 
of 10 per cent. 

In the present study, we use the redshift probabil- 
ity distributions provided by the two photometric redshift 
codes, as a quantitative tool to estimate cluster member- 
ship. Briefly, we accept galaxies as potential cluster mem- 
bers if the integrated probability for the photometric redshift 
to be within ±0.1 of the (known) cluster redshift is greater 
than a specific threshold for both of the photometric redshift 
codes. These probability thresholds (Pthresh) range from 0.1 
to 0.35, depending on the filter set available for each partic- 
ular field, and they were calibrated using our spectroscopy 
to maximise the cluster membership and, at the same time, 
to minimise contamination from interlopers. 

Calibration against our spectroscopic sample shows that 
this technique allows us to retain more than 90 per cent of 
the cluster members while rejecting slightly less than 50 per 
cent of the non-members in the spectroscopic sample. The 
efficiency of rejection for the bluest and reddest halves of the 
sample is similar to within less than 10 per cent. However, 
red non-members are slightly more efficiently rejected, and 
red members slightly more efficiently accepted than their 
blue counterparts at identical Pthresh- This is expected be- 
cause of the poorer constraints on z p hot for blue galaxies due 
to their smoother SEDs. 

A m ethod similar to that use d in our study was pro- 
posed by iBrunner fc Liibfnl (|2000h . Their application was 
based on the use of an "empirical" photometric redshift tech- 
nique. The latter is based on the use of an empirical rela- 
tion, measured for the spectroscopic sample available, be- 
tween the spectroscopic redshifts and the photometric data 
points. The use of this method, forced Brunner & Lubin to 
assume a probability distribution, which they supposed to 
be Gaussian with mean given by the estimated photomet- 
ric redshift and standard deviation defined by the estimated 
error. The error distributions are, ho wever, usually strongly 
non-Gaussian iRudnick et all 1200 ll) . so proper use of the 
probability distributions provides a better estimate of the 
real uncertainty in the photometric redshift estimates. 

When available, we use the spectroscopic information to 



determine cluster membership: spectroscopic non-members 
that are erroneously classified as cluster members by the 
photometric redshift criterion detailed above, are rejected. 
Spectroscopic members that are erroneously rejected are 
re-included into the sample before the analysis. The pho- 
tometric redshift technique we use performs well, particu- 
larly on red galaxies. As a consequence, this latter correc- 
tion based on the availability of spectroscopic information, 
does not modify significantly the results discussed below. We 
note that spec t roscop ic membership has been defined as in 
lHallidav et al] (I2004T) . The number of spectroscopic mem- 
bers for our clusters ranges from 11 to 67 (see Table 1 in 
IPoggianti et al.ll2006h . If membership is assigned using the 
photometric redshift method outlined above, the fraction of 
cluster members for which spectra are available ranges from 
a few to about 15 per cent. 

3.2 Statistical subtraction 

As an alternative method to determine cluster membership, 
we employ a 'classical' statistical subtraction technique. The 
metho d we use is similar to that adopted in lPimbblet et all 
(2002). We refer to the original paper for more details on 
the procedure which we only briefly outline here. 

The 'field' population has been determined from 
one field of the Can ada France Deep Field Survey 
jMcCracken et ahlEoOll) 3 . This corresponds to an area of 
about 0.25 deg 2 , which is much larger than the cluster area 
used for our analysis (see next section) . Both the cluster and 
the field regions are binned onto a gridded colour-magnitude 
diagram (we use a 0.3 bin in colour and a 0.5 bin in mag- 
nitude). The field region is then scaled to the same area 
as the cluster region we wish to correct, and each galaxy 
is assigned a probability to be a cluster member simply by 
counting how many galaxies lie in the colour-magnitude bin 
in the two different regions. Using a Monte Carlo method, 
the field population is then subtracted off. If this procedure 
gives a negative number of galaxies in the cluster popula- 
tion at a particular grid position, the mesh size is increased 
for that particu l ar po sition. As explained in Appendix A of 
IPimbblet e t al. (2002), this approach has the advantage of 
preserving the or iginal pr obability distribution better than 
similar methods dKodama fc B ower 200^) where the excess 
probability is distributed evenly between the neighbours of 
critical grid positions. For each cluster we run 100 Monte 
Carlo realizations of the above procedure. 

When available, we use the spectroscopic information: 
spectroscopic members and non-members are always as- 
signed a probability 1 and to be cluster members respec- 
tively. As is the case when cluster membership is assigned 
using photometric redshifts, this correction does not modify 
significantly the results presented below. 



4 THE COLOUR-MAGNITUDE RELATION 

In Fig. Qwe show the colour-magnitude diagrams for the 18 
fields for which we have high quality spectroscopy. Clusters 
are shown in order of increasing redshift. Empty circles show 

3 The catalogue has been kindly provided to us by H. McCracken 
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Figure 1. Colour-magnitude diagrams for the 18 EDisCS fields used in this study. Empty circles show objects retained by our photometric 
redshift criterion. Blue and red filled circles represent spectroscopically confirmed members with and without emission lines in their 
spectra. Thin slanted lines correspond to the 1, 3 and 5<r detection limits in the V-band. The solid thick line in each panel, represents 
the best fit relation measured using the bi-weight estimator and assuming a fixed slope of —0.09. Dashed lines correspond to ± 0.3 mag 
from the best fit line. Crosses show the location of galaxy models with two different SF histories. The solid vertical line in each panel 
shows the apparent magnitude which translates to My = —18.2 when evolved passively to z = 0. (See text for details). 
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Figure 1. — continued 



objects for which our photometric redshift criterion gives a 
high probability of cluster membership (see Sec. 13.11 . Red 
and blue filled circles represent spectroscopically confirmed 
members lacking/showing any emission line in their spec- 
tra respectively. We note that the typical detection limit 
for t he [OII13727 line is low in our spectra - approximately 
2A jPoggianti et aUl2006ft . Thin slanted lines correspond 
to the 1, 3 and 5cr detection limits in the V-band, while 
the solid thick line in each panel represents the best fit to 
the colour-magnitude relation measured using a fixed slope 
of —0.09. The fit has been computed usin g the bi-weight 
estimator feeers. Flvnn fc Gebhardl ll99Ch on the objects 
without emission lines in their spectra (red filled circles). 
Dashed lines correspond to ± 0.3 mag from the best fit line. 

Crosses in Fig. show the location of galaxy mod- 
els with two different star formation (SF) histories: a sin- 
gle burst at z = 3 (dark red) and an exponentially de- 
clining SF starting at z = 3 (green) with a characteris- 
tic time-scale (r) of 1 Gyr. Both mod els were calculated 
using the population synthesis code bv lBruzual fc Charlod 
(2003) with a Chabrier initial mass function. For each SF 
history, three different metallicities are shown: 0.02, 0.008 
and 0.004, from brighter to fainter objects. The relation be- 
tween metallicity and luminosity in these models has been 
calibrated by requiring that they reproduce the observed 
CMR in Coma (see Sec. [SJ . This calibration has been found 
a posteriori to be in good agreement with the metallicity- 
luminosit y relation derived from spectral indices of Coma 
galaxies fPoggianti et al.ll2001rl . The solid vertical line in 



each panel of Fig. Q shows the apparent magnitude which 
translates to Mv = —18.2 when evolved passively to z — 
using the single burst model. This magnitude limit was cho- 
sen so that all brighter objects are above the 5a detection 
limit in the V-band in all of our clusters. 

Note that the FORS2 field is 6'8 x 6!8 with a pixel size 
of O'.^O and, after dithering, the field of view with the maxi- 
mum depth of exposure is approximately 6.5' x 6.5'. For our 
infrared data, however, taking into account dithering and 
overlapping exposures, the field is 6.0' x 4.2' for our inter- 
mediate redshift cluster candidates and 5.4' x 4.2' for our 
high redshift cluster candidates. As explained in Sec. 13.11 
galaxies likely to be cluster members are selected as hav- 
ing a probability within ±0.1 of the cluster redshift above 
some threshold. Where there is no IR data, however, the 
redshift probability distributions are broader, and so clus- 
ter likelihoods are correspondingly lower and fewer galaxies 
meet the adopted criteria. While we have re-calibrated the 
photo-z cutoffs in the regions where no IR data are available 
so to use broader cuts in these regions, we have noted that 
some 'edge effects' remain. For these reasons, in this study 
we will only use the region for which we have both optical 
and infrared data for each of our fields. 

The open circles shown in Fig. Q correspond to objects 
within the fixed maximum physical radius centred on the 
BCG and included in the SOFI field. This physical radius 
turns out to be ~ 0.74 Mpc. For two of the fields shown 
in Fig. [T] (cll227.9-1138 and clll03.7-1245a), the BCG lies 
close to the edge of the chip (see Fig. 1 in Poggianti et al. 
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Figure 2. Evolution of the zero-point of the colour-magnitude 
relation as a function of redshift. Filled circles with error bars 
represent the values measured by fitting the observed relation. 
The solid and dashed lines show the predictions of a single burst 
model with formation redshift 3 and 2 respectively. In both cases, 
the zero-point has been measured at the apparent magnitude that 
corresponds to My = —20. when evolved to z = 0. 



2006 and Fig. 6 in White et al. 2005). In these cases, the 
open circles show all the objects with high probability of 
cluster membership within the whole region for which we 
have infrared data. 

Fig.0sh.ows several interesting results. The single burst 
model seems to provide a good fit to the observed red se- 
quence over all the redshift range sampled by the clusters 
under investigation. This confirms that the location of the 
CMR observed in distant clusters, requires high redshifts of 
formation, and that the slope is consistent with a correlation 
between galaxy metal content and luminosity. In Fig. [5] wc 
show the evolution of the zero-point of the colour-magnitude 
relation as predicted by adopting the single burst model with 
formation redshift 3 (solid line) , and the values measured by 
fitting the observed relation for all the clusters in our sample 
(filled symbols with error bars). Both for the observational 
data and for the model predictions, the zero-point has been 
measured at the apparent magnitude that corresponds to 
My — —20 when evolved to z = 0. As mentioned before, the 
fit has been obtained using only the spectroscopically con- 
firmed members without emission lines in their spectra (red 
filled circles in Fig. and assuming a fixed slope of —0.09. 
Error bars have been estimated via bootstrap resampling of 
the galaxies used to compute the fit. Overall, the measured 
values follow nicely the model line, although some devia- 
tions are visible where the models are systematically redder 
than the best fit relations (see also Fig0. These could in- 
dicate a lower formation redshift for red-sequence galaxies 
in these clusters. For a single burst model with formation 
redshift 2, the predicted zero-point is lower at all redshifts, 
as indicated by the dashed line in Fig. 01 Within the errors, 
however, both single burst models provide a relatively good 
fit over all the redshift range under investigation. 

Anot her interesti n g resu lt shown in Fig. Q was already 
noted in IWhite et alJ (|2005): our sample shows large vari- 
ations in the relative proportions of red and blue galaxies. 
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Figure 3. Fraction of galaxies bluer than 0.3 mag from the best fit 
colour magnitude relation as a function of redshift. Red and blue 
symbols are for galaxy clusters with velocity dispersion larger and 
smaller than 600 km s -1 respectively. Errors have been estimated 
using Poisson statistics. The solid and dashed line show linear fits 
to the data obtained by selecting cluster members using their pho- 
tometric redshifts and statistical subtraction respectively. Points 
are shown for the former selection criterion. (See text for details). 



Some clusters exhibit a strong red sequence with rather few 
blue galaxies (e.g. cll232. 5-1250, clll38. 2-1133) while others 
show many blue galaxies but relatively few passively evolv- 
ing systems (e.g. cll040. 7-1155, cll227.9-1138). So although 
the locus of red galaxies is well described by a uniformly 
old stellar population, the "morphology" of the colour- 
magnitude diagrams is quite varied. This must be related, 
at some level, to the dynamical and accretion histories of 
the clusters. 

In Fig.[3]we show the fraction of blue galaxies for the 18 
EDisCS fields shown in Fig.Qas a function of redshift. Red 
and blue symbols are for clusters with velocity dispersion 
larger and smaller than 600 km s -1 respectively. The blue 
fractions shown in Fig. |3] have been computed by count- 
ing the galaxies bluer than 0.3 from the best fit colour- 
magnitude relation, and brighter than the passive evolution 
corrected limit that corresponds to —18.2 in the rest-frame 
V-band at z = 0. All the objects which have a high proba- 
bility of cluster-membership and within the same areas used 
for Fig. Q are used. The solid line shows a linear fit to the 
data, while the dashed line shows a linear fit to data ob- 
tained using a statistical subtraction technique to determine 
cluster membership. Errors have been estimated using Pois- 
son statistics. Our definition of b lue fraction differs from th e 
original definition introduced bv lButcher fc Oemlerl lll984l) . 
It is, however, interesting to note that res ults shown in Fig. El 
exhibit - as found for the first time by iButcher fc Oemlerl 
( 1984) - an increase in the fraction of blue galaxies with in- 
creasing redshift. The trend is present, although weaker, also 
when a statistical field subtraction instead of photometric 
redshift information is used to determine cluster member- 
ship. A Spearman's rank correlation test gives a coefficient 
of 0.50 with a significance level of 0.036 in the case mem- 
bership is defined using statistical subtraction, while in the 
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case membership is defined using photometric redshifts, the 
coefficient is 0.69 with a significance level of 0.001. It should 
be noted, however, that the error bars are quite large and 
that there are large cluster-to-cluster variations at any given 
redshift. 

We note that the constraints on the photometric red- 
shift for blue galaxies are usually worse than those for galax- 
ies of the same magnitude but with a redder colour, because 
of their smoother spectral energy distribution. The statis- 
tical subtraction technique is also more uncertain for blue 
galaxies because the field population that is used to perform 
the subtraction has a colour distribution peaked towards 
blue colours. In the following sections, we will concentrate 
on the distribution of galaxies along the red-sequence, where 
both the photometric redshift and the statistical subtraction 
method are expected to perform better. 



5 THE RED-SEQUENCE GALAXY 
DISTRIBUTION 



In iDe Lucia et al.1 (I2004I) . we analysed the distribution of 
galaxies along the colour-magnitude relation for four of the 
highest redshift clusters in the EDisCS sample. As men- 
tioned in Sec.0 our analysis pointed out a significant deficit 
of faint red galaxies compared to the nearby Coma cluster. 
We interpreted this deficit as evidence for a large number of 
galaxies having moved onto the red-sequence relatively re- 
cently, possibly as a consequence of the suppression of their 
star formation by the dense cluster environment. If the sce- 
nario we envisaged is correct, then we should be able to 
see some evolution in the relative number of 'faint' and 'lu- 
minous' red galaxies as a function of redshift within our 
EDisCS sample. 

In Fig. [I] we show the number of galaxies along the 
red-sequence obtained by averaging the distributions of in- 
dividual clusters in three redshift bins. Black histograms are 
obtained by selecting all the galaxies for which our photo- 
metric redshift criterion gives a high probability of cluster 
membership (Sec. l3.il . Red histograms are obtained by se- 
lecting cluster members using a purely statistical subtrac- 
tion (Sec. YS.2) . In the latter case, we have averaged over 
100 Monte Carlo realizations of the statistical subtraction 
for each cluster. For the histograms shown in Fig. |1J all the 
'members' within ~ 0.5 x R200 and within ±0.3 mag of 
the best fit relation are used. In the present analysis, we 
are neglecting the clusters cll227.9-1138 and clll03.7-1245a 
for which the BCG lies at the edge of the chip. Furthermore, 
w e also neglect t he clu ster cllll9. 3-1129 which, as explained 
in lWhite et alJ (120051) . shows a very weak concentration of 
galaxies close to the selected BCG and has a small value of 
R2oo- As a consequence, there are just a handful of galaxies 
on the red sequence within the fraction of virial radius we 
have adopted. In addition, we do not have infrared data for 
this cluster. 

Clusters have been combined in three redshift bins (we 
end up with 5 clusters in each bin) correcting colours and 
magnitudes to the central redshift of the corresponding bin. 
Corrections are based on the single burst model shown in 
Fig. Q (results do not appreciably change using a single burst 
model with formation redshift 2 instead of 3). The scale on 
the top of each panel in Fig.0]shows the rest-frame V-band 
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Figure 4. Number of galaxies along the red— sequence. His- 
tograms have been obtained by stacking individual clusters in 
three redshift bins. Black and red histograms correspond to clus- 
ter membership based on photometric redshift and on statisti- 
cal subtraction respectively. The scale on the top of each panel 
shows the rest— frame V-band magnitude that corresponds to the 
I-band magnitude and has been passively evolved to z = 0. Verti- 
cal dashed lines show our magnitude limit and the edge between 
'luminous' and 'faint' galaxies. A small offset has been added 
between two histograms in each panel for clarity. For each his- 
togram, we give the luminous— to-faint ratio as defined in the 
text. 



magnitude corresponding to the observed I-band magnitude 
after correcting for passive evolution between the redshift of 

the bin a nd z = 0. 

As in lDe Lucia et all i2004f) . we compute for each red- 
shift bin a 'luminous-to-faint' ratio. We classify as 'lumi- 
nous' all galaxies brighter than My = —20.0 and as 'faint' 
those galaxies that are fainter than this magnitude and 
brighter than My — —18.2. As mentioned in the previous 
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Figure 5. Luminous— to-faint ratio as a function of redshift. Cir- 
cles and triangles are used in the case where all galaxies within 
~ 0.5 X R200 or within ~ 0.74 Mpc from the BCG are selected. 
Black and red symbols correspond to cluster membership based 
on photometric redshift and on statistical subtraction. Symbols 
corresponding to the same redshift have been displaced for clarity. 



section, the faint limit has been chosen because it corre- 
sponds to the limiting magnitude in the I-band for which 
all selected objects are above the 5a detection limit in the 
V-band. As for the choice of the magnitude corresponding 
to the edge between faint and luminous galaxies, we use —20 
because, at our highest redshift, it approximately equally di- 
vides the range of cluster galaxy magnitudes covered down 
to —18.2. Both limits correspond to values after passive evo- 
lution to z = and are indicated by vertical dashed lines in 
Fig. El 

The values of the luminous-to-faint ratios computed 
for the histograms shown in Fig. 21 are listed in each panel, 
along with the estimated errors. Fig. |5] shows these values 
as a function of redshift. The error bars have been esti- 
mated assuming Poisson statistics and, in the case where 
cluster members are selected using a statistical subtraction, 
they include the error contribution from the background 
field (this is however negligible given the large area used 
for the subtraction). Circles in Fig. |5] correspond to the his- 
tograms plotted in Fig. [I] where all the cluster members 
within ~ 0.5 x R200 have been used. Triangles are used in 
the case where all the cluster members within a fixed phys- 
ical distance (~ 0.74 Mpc) from the BCG are retained. Red 
and black symbols correspond to membership criteria based 
on statistical subtraction and on photometric redshifts re- 
spectively. 

The error bars shown in Fig.[H]are large, and some small 
differences arise from the use of different criteria for clus- 
ter membership and from different choices about the area 
used for the analysis. Overall, however, independently of 
the method employed and the area used, the data indicate 
a decrease of the luminous-to-faint ratio with decreasing 
redshift. Faint red galaxies become increasingly important 
with decreasing redshift or, in other words, the faint end 
of the colour-magnitude relation becomes increasingly pop- 
ulated with decreasing redshift. As noted in our previous 



paper, this finding is inconsistent with a formation scenario 
in which all red galaxies in clusters today evolved passively 
after a synchronous short duration event at z > 2 — 3, and 
suggests that present day passive galaxies follow different 
evolutionary paths, depending on their luminosity. 

It is now interesting to ask if this evolution in the 
luminous-to-faint ratio depends on cluster properties, for 
example mass or velocity dispersion. In Fig.|S|we again show 
the distribution of galaxies along the red-sequence. This 
time, we have combined the clusters in two redshift bins 
and, in each redshift bin, we have split the clusters accord- 
ing to their velocity dispersions. Red and blue histograms 
are for clusters with velocity dispersions larger and smaller 
than 600 kms -1 respectively. Black histograms are obtained 
by stacking all the clusters in each redshift bin. Left panels 
are for the case where membership is based on photomet- 
ric redshifts, while for the histograms shown on the right, 
membership is based on a purely statistical subtraction. All 
cluster members within ~ 0.5 x R200 from the BCG are used. 
In the lower redshift bin, we have 5 clusters in each bin of 
velocity dispersion while in the higher redshift bin we have 
2 and 3 clusters in the larger and smaller velocity disper- 
sion bin respectively. The behaviour shown in Fig. |S| is not 
significantly different if all galaxies within a fixed physical 
radius from the BCG are used (see also Fig. 01. 

The values listed in Fig. |S| show the same trend of an 
increasing luminous-to-faint ratio as a function of redshift 
and also hint at a dependence on cluster velocity dispersion. 
Clusters with large velocity dispersion seem to have a larger 
fraction of luminous galaxies with respect to the systems 
with smaller velocity dispersion. For the highest redshift bin, 
the difference goes in the same direction but is not statis- 
tically significant. The number statistics are, however, poor 
and the error bars are large so that it is difficult to draw 
any definitive conclusions. We find, however, similar results 
if we split the cluster s on the basis of a richness estimate 
similar to that used in lWhite et alJ (l2005h . i.e. based on the 
number of red-sequence galaxies. 

In a recent study, Ta naka et alJ {2005) have investi- 
gated the build up of the colour-magnitude relation using 
deep panoramic imaging of two clusters at z = 0.83 and 
2 = 0.55 respectively. Using photometric redshifts and sta- 
tistical subtraction, and using nearest-neighbour density to 
characterise the environment, these authors conclude that 
build-up of the colour-magnitude relation is "delayed" in 
lower density environments. This is in apparent contradic- 
tion with our findings, although a direct comparison is not 
straightforward as we use the cluster velocity dispersion and 
not local density t o differentiate e nviron ments. In addition, 
the conclusions of iTanaka et alJ {2005) are based only on 
two clusters and these authors argue that their intermedi- 
ate redshift cluster is "peculiar" . Finally, one should keep in 
mind that cluster-to-cluster variations are rather large (see 
Fig. 0. Further studies are therefore needed to confirm or 
disprove the apparent trends. 



6 THE RED-SEQUENCE GALAXY 

DISTRIBUTION IN NEARBY CLUSTERS 

In the previous section, we have analysed the evolution of 
the luminous-to-faint ratio over the redshift range sampled 
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Figure 6. As in Fig. [I]for two redshift bins. Red and blue histograms are for clusters with velocity dispersion larger and smaller than 
600kms — 1 respectively. Black histograms are obtained by stacking all the clusters in each redshift bin. Left and right panels correspond 
to membership based on photometric redshifts and on statistical subtraction respectively. All cluster members within ~ 0.5 X R200 are 
used. Histograms in each panel have been slightly offset for clarity. The luminous— to-faint ratio (with the corresponding errors) are listed 
in each panel. 



by our EDisCS clusters. We want now to set the zero-point 
for this evolution by studying the distribution of galaxies 
along the red-sequence in nearby galaxy clusters. In order 
to carry out a comparative study with the clusters in the 
EDisCS sample, we need relatively deep photometry (down 
to —18.2 in the rest-frame V-band) sampling the rest-frame 
U and V bands, and with good spatial coverage (at least 
~ 0.5 x R,2uo)- These conditions are met by only a few nearby 
clusters. In the following we will describe in more detail 
the data and the analysis performed on our low redshift 
comparison sample, which is constituted by Coma, and a 
sample of clusters selected from the Sloan Digital Sky Survey 
(SDSS). 

To date, Coma (A1656) remains the only rich cluster 
in the nearby Universe (2 = 0.023) with a high precision 
near-UV colour-magnitude relation, determined using hun- 



dreds of spectroscopically confirmed members. In this work, 
we have used magnitudes and colours in a 25'/ 2 diameter 
aperture from iTerlevich. Caldwell fc Bowerl l|200lft . At the 
redshift of Coma, this corresponds to a physical size of 11.71 
kpc, quite closely approximating our ~ 11-15 kpc aperture 
from z ~ 0.4-0.8 (see Sec.^J. Observed magnitudes are con- 
verted to absolute magnitudes using the distance modulus of 
Coma computed using its redshift and the adopted cosmol- 
ogy (35.00). Observed colours are co nverted to rest- frame 
colours using tabulated K-corrections llPoggiantll997l) . The 
colour-magnitu de relation of the Com a cluster, based on the 
catalogue bv ITerlevich et al] ll200ll) . is plotted in the top 
panel of Fig. [7| The solid thick line shows the best fit re- 
lation determined using the bi-weight estimator, as for the 
EDisCS clusters. Dashed lines correspond to ±0.3 mag from 
the best fit line. Crosses show predictions from the same two 
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Figure 7. Colour-magnitude relation (top panel) and distribu- 
tion of galaxies along the red-sequence (bottom panel) for the 
Coma cluster. The solid line in the top panel shows the best fit 
relation to the red— sequence. Crosses correspond to the same SF 
models shown in Fig. HI 



models shown in Fig. Q We recall that the relation between 
metallicity and luminosity in these models has been cali- 
brated by requiring that they reproduce the observed CMR 
in Coma. 

The bottom panel of Fig.Qshows the distribution along 
the red-sequence for this cluster. Membership information 
has been obtained using a redshift catalogue kindly pro- 
vided by M atthew Colless and a p rocedure similar to that 
employed in lMobasher et all (|2003f) ■ Briefly, for each magni- 
tude bin, we count how many objects have a measured spec- 
troscopic redshift (N z ), and how many are spectroscopically 
confirmed members (N c ). We assume then that the spec- 
troscopic sample is 'representative', i.e. that the fraction of 
galaxies that are cluster members is the same in the spec- 
troscopic sample (that is incomplete) as in the photometric 
sample (that is complete). Cluster membership can then be 
obtained as the ratio between the two numbers computed 
before (N c /N z ). The counts shown in the bottom panel of 
Fig. |7] are obtained correcting the raw distribution by this 
membership factor. The luminous-to-faint ratio we measure 
for Coma is 0.32 ± 0.06. 4 
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Figure 8. Stacked colour— magnitude relation for galaxy clusters 
selected from the C4 catalogue with velocity dispersion smaller 
(top panel) and larger than 600kms — 1 . The solid line in each 
panel shows the colour— magnitude relation predicted by a single 
burst model with formation redshift 3. Dashed lines are offset 
from the solid line by ±0.1. The vertical dashed lines correspond 
to M v = -20.18 and M v = -18.25 respectively. (See text for 
details). 



We have complemented our low redshift comparison 
sample with a sample of clusters selected from the SDSS. 
The basis fo r the cluster s ample used here is the C4 cluster 
catalogue bv lMiller e t al. (2005). This catalogue is available 
for the SDSS Data Release 3 and is based on the spectro- 
scopic sample. The cluster detection algorithm employed for 
the construction of the C4 catalogue, is essentially based on 
an identification in a three-dimensional space (position, red- 
shift, and colour). We refer to the original paper for more de- 
tails. Based on the cluster redshifts and velocity dispersions 
given in the C4 catalogue, we have re-identified the BCG 
and measured the velocity dispersion at the virial radius for 
each cluster. Details about the procedure are described in 
von der Linden et al. (in preparation). For the purposes of 
this analysis, we have used the 76 clusters below z = 0.045 



4 In lDe Lucia et all 12004F) we erroneously used a distance mod- 
ulus equal to 35.16 (instead of 35.00 used here) corresponding to 



Ho = 65kms~ 1 Mpc - 1 . However, a small bug in the code we 
used produced a value of the luminous— to-faint ratio that is not 
much off that measured in the present study (0.34 ± 0.06). 
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so to assure completeness down to the magnitude limit used 
for our analysis. 

A direct comparison with the analysis of the red- 
sequence galaxies distribution performed above for the 
EDisCS clusters and for Coma is not simple, and requires a 
number of steps that we describe in the following. The first 
difficulty comes from the fact that, for the SDSS clusters, 
we have AB Petrosian magnitudes, while for the EDisCS 
clusters and for Coma we have used aperture magnitudes. 
In order to have an estimate of the correction necessary to 
convert Petrosian magnitudes into aperture magnitudes, we 
have compared our I-band aperture mag n itude s to the 'to- 
tal' magnitudes we used in IWhite et all {2005). We recall 
that an approximate total I-magnitude for each galaxy was 
estimated by adding to the Kron magnitude the correction 
appropriate for a point source measured within an aperture 
equal to the galaxy's Kron aperture (we refer to the original 
paper for details). Using the median value of this correction, 
and considering that for elliptical like objects Petrosian mag- 
nitudes take into account about 80 per cent of the light, we 
estimate that our aperture magnitudes can be converted into 
Petrosian magnitudes by a shift that varies between —0.05 
and —0.18 from fainter to brighter galaxies. 



In practice, we have converted the limits used before 
(M v = -18.2 and M v = -20) into "Petrosian limits" 
(Mv = -18.25 and M v = -20.18). The conversion from 
apparent magnitudes to absolute magnitudes in the V- 
band has been perfor med using the routine KCORRECT by 
iBlanton et alj J2003f) . For each cluster we have then con- 
structed photometric catalogues by taking all the objects 
that are classified as galaxies by the SDSS pipeline and that 
reside within ~ 0.5 x R200 from the BCG. A corresponding 
field catalogue has been constructed by using the whole DR3 
and for each cluster we have performed 50 Monte Carlo re- 
alizations of the statistical subtraction procedure described 
in Sec. |3~21 

Fig. |H] shows, for one realization of the statistical sub- 
traction, the stacked colour-magnitude diagrams for the two 
velocity dispersion bins. The solid line in each panel shows 
the relation predicted by a single burst model with forma- 
tion redshift 3, while dashed lines are offset from this line 
by ±0.1 mag. The vertical dashed lines show the limits used 
to define the luminous and faint population. We note that 
the magnitudes and colours plotted in Fig. |S| are given in 
the AB system. We have therefore used AB magnitudes and 
colours for the single burst model. The luminous-to-faint 
ratio computed from these stacked colour-magnitude rela- 
tions are listed in each panel and are 0.39 ±0.04 for the high 
velocity dispersion bin, and 0.50 ± 0.03 for the low veloc- 
ity dispersion bin. The value measured for the clusters with 
velocity dispersion larger than 600 kms -1 appears then com- 
patible, within the errors, with that measured for the Coma 
cluster, while the value measured for clusters with lower ve- 
locity dispersion is significantly higher. We note also that the 
trend found for the SDSS clusters is the opposite of what 
we have found for our high redshift sample. 



7 THE BUILDUP OF THE 

COLOUR-MAGNITUDE RELATION 

The results of our analysis are summarised in Fig. [5] Filled 
black circles show the values of the luminous-to-faint ratio 
for EDisCS clusters in three redshift bins. The values shown 
in Fig.0 obtained for different choices of cluster membership 
criterion and area, have been averaged together. The error 
bars corresponding to these points have also been averaged, 
rather than combined in quadrature, in order to give a 'con- 
servative' measure for the uncertainties. The green diamond 
shows the corresponding value for Coma. The orange and 
cyan triangles show the value measured for clusters selected 
from the SDSS with velocity dispersion larger and smaller 
than 600 kms -1 respectively. 

In iDe Lucia et al] {20Q4 1 ), we interpreted the deficit of 
faint galaxies found in the high redshift EDisCS clusters, as 
evidence that a large fraction of present-day passive faint 
galaxies must have moved onto the colour-magnitude rela- 
tion at redshift lower than 0.8. We argued that the popu- 
lation of blue galaxies observed in distant galaxy clusters 
provide the logical progenitors of faint red galaxies at z = 0. 
It is therefore interesting to ask if the measured evolution 
in the luminous-to-faint ratio can be reproduced by simple 
evolution of the combined blue and red galaxies that pop- 
ulate the colour-magnitude diagrams of our high redshift 
clusters. In order to have a handle on th is question, we have 
used t he population synthesis model bv lBruzual fc Charlod 
( 2003) to construct different exponentially declining star for- 
mation histories with r = 1, 2, 3, and 7 Gyr and a redshift of 
formation 3. The same metallicities and normalisations used 
for the models shown in Fig. Q have been adopted here. We 
have then started from the distribution of galaxies on the 
colour-magnitude diagram of the clusters in the highest red- 
shift bin shown in Fig.|U]- for simplicity, we have used mem- 
bership based on photometric redshift and a fixed physical 
distance from the BCG. For each galaxy bluer than 0.3 mag 
than the best fit relation, we have determined the "closest 
model" in colour-magnitude space among those listed above. 
Each galaxy is then evolved to the next redshift bin by the 
amount predicted by the best fit model with the correspond- 
ing star formation history truncated at the redshift of the 
observation of the cluster (truncation model) or 1 Gyr later 
(delayed model). For the red-sequence galaxies (those in the 
stripe used to compute the luminous-to-faint ratio) , we have 
simply used the single burst model described in Sec. 2] Us- 
ing the same method, we use the distribution of galaxies in 
the colour-magnitude diagram of all higher redshift clusters 
to predict the evolution to z — 0.45 and to z — 0.025. 

Model predictions are shown as red filled (truncation 
model) and open blue (delayed model) squares connected 
by arrows in Fig. [5] Interestingly, both models predict an 
amount of evolution from 2 = 0.75 to z ~ that is in 
nice agreement with that measured using the highest red- 
shift EDisCS clusters and the Coma cluster. The predicted 
amount of evolution between z ~ 0.4 and z ~ is instead 
too large to reproduce the luminous-to-faint ratios mea- 
sured for the SDSS clusters. These appear to be compatible 
with the luminous-to-faint ratio measured for the EDisCS 
clusters at z ~ 0.4. In the truncation model, galaxies move 
onto the red-sequence very rapidly so that the predicted 
luminous-to-faint ratio lies below the measured value in the 
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Figure 9. Luminous— to-faint ratio as a function of redshift. 
Filled circles show the average of the values shown in Fig.|3 The 
green diamond shows the values determined for Coma. The trian- 
gles show results for clusters selected from the SDSS with velocity 
dispersion larger (orange) and smaller (cyan) than 600kms — 1 . 
The arrows indicate the evolution of the luminous-to-faint ratio 
obtained using the models described in the text. 



intermediate redshift bins. In the delayed model, galaxies 
stay blue longer so that the predicted evolution is closer to 
the observed trend at all redshifts sampled by the EDisCS 
clusters. 

In Sec. |K| we have investigated how the distribution of 
galaxies on the red-sequence depends on the cluster ve- 
locity dispersion. The corresponding luminous-to-faint ra- 
tios (again averaged for different cluster membership crite- 
ria and different choices for the area) are shown in both 
panels of Fig. ED Filled red and open blue circles are used 
for clusters with velocity dispersion larger and smaller than 
600 km s -1 respectively. The arrows connected by squares 
show the evolution predicted by the truncation (top panel) 
and delay (bottom panel) models described above. Triangles 
refer to the SDSS clusters as in Fig.UD Both models predict 
a luminous-to-faint ratio that is close to observed value for 
low velocity dispersion clusters at redshift z ~ 0.5. For clus- 
ters with velocity dispersion larger than 600kms _1 , both 
models instead predict a lower value than that measured for 
intermediate redshift EDisCS clusters. The values predicted 
for clusters with larger velocity dispersion are in both mod- 
els and down to redshift zero, larger than the corresponding 
values for low velocity dispersion clusters. This appears in 
agreement with the trend found within our EDisCS, but in 
contradiction with that found for the SDSS clusters. In ad- 
dition, as shown already in Fig. [5] the amount of evolution 
predicted between z ~ 0.5 and z ~ is not compatible with 
the high luminous-to-faint ratios measured for the SDSS 
clusters. 

We note that Table 4 in IPoggianti et alJ ll2006l) indi- 
cates that a cluster with velocity dispersion ~ 600 km s" 1 
at z ~ 0.6 evolves into a system with velocity dispersion 
~ 700 km s" 1 at z ~ 0. Therefore, it would be more correct 




Figure 10. Luminous— to-faint ratio as a function of redshift for 
clusters split according to their velocity dispersions. Filled red 
and open blue circles are used for EDisCS clusters with velocity 
dispersions larger and smaller than 600 km s -1 respectively. Tri- 
angles refer to the SDSS clusters as in Fig. [5] Squares connected 
by arrows show prediction from the truncation (top panel) and 
delay (bottom panel) described in the text. 



to compare the values predicted from our models at z ~ 
with those obtained for the SDSS by using a cut correspond- 
ing to ~ 700 km s . In this case, the measured luminous- 
to-faint ratios are 0.47 ± 0.03 and 0.44 ± 0.06 for the low 
and high velocity dispersion clusters respectively. The sam- 
ple used in this study, however, only contains 4 cluster with 
velocity dispersion larger than ~ 700 km s -1 . 



8 DISCUSSION AND CONCLUSIONS 

The tight colour-magnitude relation shown by cluster ellip- 
tical galaxies has been the subject of numerous studies that, 
in the last decade, have been pushed to higher and higher 
redshifts. An interesting and still controversial result of re- 
cent years is the claim of a 'deficit' of faint red galaxies a t 
redshift z ~ 1. In our previous work jDe Lucia et al]l2004l) . 
we analysed the distribution of galaxies on the red-sequence 
for four high redshift EDisCS clusters, and we compared 
such a distribution to that measured for the nearby Coma 
cluster. Although with a low significance level, we found 
that clusters at redshift z ~ 0.8 exhibit a 'deficit' of faint 
red sequence galaxies. 

A decline in the number of red sequence members at 
fai nt m agnitud es wa s first observed in clusters at z = 0.25 
bv lSmail et alJ (' 1998 ) . Evidence for a 'truncation' of the red- 
sequence was noted in an overdensity around a radio galaxy 
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at z = 1.2 by IKaiisawa et al.l (l2000f) and INakata et alJ 
(2001). The same authors, however, speculated that their 
result might be spurious because of limited area coverage 
(< 0.33 Mpc) a nd strong luminosity s egregation. (See also 
the discussion bv lKodama et aP J2004T) who obtain a similar 
result for early-type galaxies i n a sin gle deep field). 

In recent work, lAndreoP l|200fJ l has studied the red- 
sequence luminosity function for the cluster MS1054-0321. 
By comparing his result with the faint -end slope measured 
by fitting the red-sequence luminosity function of nearby 
clusters from the SDSS by iTanaka et al.l (120051) . he con- 
cluded that there is no evidence for a decreasing number 
of faint red galaxies at higher redshift. The results obtained 
by Andreon are based on a single cluster. The fitting proce- 
dure he adopted for the luminosity function of MS1054-0321 
is not the same as that employed by Tanaka et al., and a 
comparison based on best fit parameters is plagued by the 
well known covariance between errors on M, and on a. In 
addition, we note that Tanaka et al. use an estimate of lo- 
cal density based on nearest-neighbour statistics. This com- 
plicates the comparison with MS1054-0321, a cluster with 
large velocity dispe rsion, bright X-ray luminosity and ev i- 
dent substructures JPonahue et al]ll998l : iTran et al]ll99Sj) . 

Evolution in the distribution of red-sequence galax- 
ies is expected as a natural consequence of the recently 
established mass- dependence of ellipt ical galaxy evolution 
both in clusters and in the field j Thom^s e^alj [2QTJ5J: 
lOallazzi et al.1 Eooi Ivan der Wei et al I 1200,4 iTreu et*al] 
120051: iHolden et all 120051 : iNelan et al1l2005l) . The discus- 
sion above suggests the details are still uncertain. The er- 
ror bars are large and large cluster-to-cluster variations 
preclude definitive conclusions. Nevertheless, the build-up 
of the colour-magnitude relation is of great interest, as it 
can constrain the relative importance of star formation and 
metallicity in establishing the observed properties of ellipti- 
cal galaxies. 

In this paper, we have used galaxy clusters from our 
EDisCS sample in order to extend our previous analysis to 
a wider redshift range and to study how this effect depends 
on cluster velocity dispersion. In agreement with previous 
work, we find that bright red-sequence galaxies in high red- 
shift clusters can be described as an old, passively-evolving 
population. A single burst model with formation redshift of 
3, calibrated on the colour-magnitude relation of the nearby 
Coma cluster, provides a good fit to the red sequence ob- 
served over the full redshift range sampled by our clusters. 
This confirms earlier claims that the location of the red- 
sequence in distant clusters suggests high formation redshift, 
and that its slope is consistent with a correlation between 
metal content and luminosity. 

However, within the same EDisCS sample, we also con- 
firm our previous finding of a significant evolution in the lu- 
minosity distribution of red-sequence galaxies since z ~ 0.8. 
Combining clusters in three different redshift bins, and defin- 
ing as 'faint' all galaxies in the passive evolution corrected 
range 0.4 >L/L„ > 0.1, we find a clear decrease in the 
luminous-to-faint ratio with decreasing redshift. The error 
bars and the cluster-to-cluster variation are large, but the 
measured trend is robust against variations in the criteria 
adopted for cluster membership and in the size of the re- 
gion analysed. We have also investigated how this evolution 
depends on cluster velocity dispersion. At intermediate red- 



shift, the luminous-to-faint ratio of clusters with velocity 
dispersion larger than 600kms _1 appears to be larger than 
that measured for clusters at the same redshift but with 
lower velocity dispersion. The error bars and the cluster- 
to-cluster variations are, however, too large to draw any 
definitive conclusions regarding this point. 

Our low redshift comparison sample includes the Coma 
cluster, and a sample of clusters selected from the SDSS. For 
the Coma cluster, we find a value of the luminous-to-faint 
ratio that is significantly lower than the value obtained for 
the EDisCS clusters at z ~ 0.45. This is not the case for the 
luminous-to-faint ratios measured by stacking clusters from 
the SDSS in different velocity dispersion bins. These values 
are not significantly different from the values measured for 
the EDisCS clusters at z ~ 0.45. 

Interestingly, we find the measured amount of evolu- 
tion in the luminous-to-faint ratio from z = 0.75 to z ~ 
to be approximatively consistent with predictions of sim- 
ple models where the blue bright galaxies that populate the 
colour-magnitude diagram of high redshift clusters, have 
their star formation truncated by the hostile cluster envi- 
ronment. Clearly the model we use is extremely simplified. 
We are assuming a single redshift of formation for all galax- 
ies and guessing their star formation history simply on the 
basis of their location in the observed colour-magnitude di- 
agram. In reality, galaxies will have a certain distribution of 
formation redshifts and this, together with age, metallicity 
and dust degeneracies, will certainly complicate the mod- 
elling. In addition, we are simply assuming that the star 
formation history is truncated at the cluster redshift (or 1 
Gyr later) for all galaxies bluer than 0.3 mag from the best 
fit red-sequence. Not all these galaxies are falling into the 
cluster at the time of our observations and the time-scale 
of the star formation suppression by the cluster environ- 
ment (if there is a suppression of the star formation by the 
cluster environment) might be different than that assumed 
and/or depend on cluster or galaxy properties. Finally, we 
are neglecting further infall of galaxies between our various 
redshift bins. For all these reasons, our model results should 
be taken with caution. They simply suggest that a scenario 
in which infalling galaxies have their star formation histories 
truncated by the hostile cluster environment is in qualita- 
tive agreement with the observed build up of the colour- 
magnitude sequence. They do not yet convincingly confirm 
this scenario. 

Our results indicate that present-day passive galaxies 
follow different evolutionary paths, depending on their lu- 
minosity (or mass). This conclusion is in line with recent 
res ults from fundamental-plane and stellar population stud- 
ies llThomas et al.ll2004lvan der Wei et al.ll2005l : lTreu et alJ 
12004 IHolden et alJl2005h . More data are required to clar- 
ify if and how this depends on environment. Such studies 
will constrain the relative importance of star formation and 
metallicity in establishing the observed red-sequence, and 
thus clarify the physical mechanisms that drive the forma- 
tion and evolution of the early type galaxy population in 
clusters. 
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